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Abstract. We consider a coupled system of Hamiltonian partial differential 
equations introduced by Popowicz, which has the appearance of a two-field 
coupling between the Camassa-Holm and Degasperis-Procesi equations. The 
latter equations are both known to be integrable, and admit peaked soliton 
(peakon) solutions with discontinuous derivatives at the peaks. A combination 
of a reciprocal transformation with Painleve analysis provides strong evidence 
that the Popowicz system is non-integrable. Nevertheless, we are able to con- 
struct exact travelling wave solutions in terms of an elliptic integral, together 
with a degenerate travelling wave corresponding to a single peakon. We also 
describe the dynamics of JV-peakon solutions, which is given in terms of a 
Hamiltonian system on a phase space of dimension 3N. 

1. Introduction. The members of a one-parameter family of partial differential 
equations, namely 

mt + Unix + bUxTn = 0, m = u — u^x (1) 

with parameter &, have been studied recently. The case 6 = 2 is the Camassa-Holm 
equation [1], while 6 = 3 is the Degasperis-Procesi equation [•']], and it is known 
that (with the possible exception of 6 = 0) these are the only integrable cases 
[13], while all of of these equations (apart from h = —1) arise as a shallow water 
approximation to the Euler equations [6] . All of the equations have at least one 
Hamiltonian structure [11], this being given by 

mt=B^-^^ B=-b^m'-'^'dxm'/'Gm'/'dxm'-'/\ (2) 
dm 

with G = {dx — dx)^^ and the Hamiltonian H = {b — 1)^^ j mdx for 6 ^ 0, 1 (and 
the latter special cases admit a similar expression). 

One of the most interesting features of these equations is that their soliton solu- 
tions are not smooth, but rather the field u has a discontinuous derivative at one or 
more peaks (hence the name peakons), while the corresponding field m is measure 
valued. More precisely for the single peakon the solution has the form 

u = c exp(a; — ct — Xq), with m = 2c 5{x ~ ct ~ xg) 
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(with a;o being an arbitrary constant), while iV-peakon sohitions are given by 

N N 

'^^ = ^Pjii) ^M~\x~<ljit)\), ■m = 2^pj{t)S{x - qj{t)), (3) 

where the amplitudes (t) and peak positions qj(t) satisfy a Hamiltonian dynamical 
system for any b. For 6 = 2 the qj and pj arc canonically conjugate position 
and momentum variables for an integrable geodesic flow with the co-metric g^'' ~ 
exp{ — \qj — qk\) [I], and for 6 = 3 the peakon motion is again integrable, being 
described by Hamilton's equations for a different Poisson structure [4, 5], but for 
arbitrary b the 7V-peakon dynamics is unlikely to be integrable in general [10]. 

There is currently much interest in generalisations of the Camassa-Holm equa- 
tions and its relatives. Qiao has found an integrable equation of this type with 
cubic nonlinearity [17], and another such example was discovered very recently by 
Vladimir Novikov [ J 4]; one of the authors spoke at the AIMS meeting in Arlington 
on this topic (for more details see [12]). An important challenge is to understand 
the solutions of coupled equations with two or components, and in higher dimen- 
sions. For example the EPDiff equation can be used to describe fluids in two or 
more spatial dimensions, as well as appearing in computational anatomy [!)]. Chen 
et al. found an integrable two-component analogue of the Camassa-Holm equation 
[2] , which also admits multi-peakon solutions [7] . Popowicz has constructed another 
two-component Camassa-Holm equation using supersymmetry algebra [l'<\. 

The purpose of this short note is to summarise some preliminary results that we 
have obtained on the two-component system given by 

mt + mx{2u + v) + 3m(2ux + Vx) = 0, 
nt + n42u + v) + 2n{2ux + vj) = 0, (4) 
m^u-Uxx, n^v-Vxx- 

This system can be considered as a coupling between the Camassa-Holm equation 
and the Degasperis-Procesi equation (corresponding to (1) for 6 = 2,3 respectively); 
it reduces to the former when u — 0, and to the latter when w = 0. The system (4) 
was obtained by Popowicz by taking a Dirac reduction of a three-field local Hamil- 
tonian operator [16]. By construction, this system has a (nonlocal) Hamiltonian 
structure, and due to the existence of conservation laws it was conjectured that it 
should be integrable (although no second Hamiltonian structure was found). 

After reviewing the Hamiltonian structure for it in the next section, in section 3 
we perform a reciprocal transformation on the system (a nonlocal change of inde- 
pendent variables) which transforms it to a third order partial differential equation 
for a single scalar field. By applying Painlevc analysis of the singularities in solu- 
tions of the reciprocally transformed system, we find the presence of logarithmic 
branching, which is a strong indicator of non-integrability. Nevertheless, in section 
4 we find that the system has exact travelling wave solutions given by an elliptic 
integral, as well as a degenerate travelling wave which is a peakon. In section 5 we 
present formulae for iV-peakon solutions of (4), which are governed by Hamiltonian 
dynamics on a 37V-dimensional phase space. The final section is devoted to some 
conclusions. 
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2. Hamiltonian and Poisson structure. Popowicz constructed the system (4) 
from the Hamiltonian operator 

9m2/3a,mi/3G'TOi/3a^^2/3 6m2/39,mi/3Gni/2a^„i/2 



(5) 



where G = {dx — d^) ^- With the Hamiltonian 

Ho = J {171 + n) dx, (6) 



the system can be written as 



nt 



Sm 

= Z\ |={m,ffo}. 

Mo 

(5n 



For Z as in (5), the Poisson bracket between two functionals A, B is given by the 
standard formula 

{A,B} = [(^ ^ W ( #) ) dz, 

J ^ ^ \ 5n{z) j 

which is equivalent to specifying the local Poisson brackets between the fields m 
and n as 

{m(x),m(j/)} = mx{x)mx{y)G{x - y) 

+3{m{x)mx{y) - mx{x)m{y))G' {x - y) 
—9m{x)m{y)G" [x ~ y), 
{■m{x),n{y)} = mx{x)nx{y)G{x - y) 

+ {3m{x)nx{y) - 2mx{x)n{y))G' {x - y) (7) 
—6m{x)n{y)G" {x — y), 
{n{x),n{y)} = nx{x)nx{y)G{x - y) 

+2{n{x)nx{y) - nx{x)n{y))G' {x - y) 
~4:n{x)n{y)G"{x~y), 

where 

Gix) = isgn(x) (l - e-l-l) (8) 

is the Green's function of the operator G. This G satisfies the functional equation 

G'{a) (g(/3) + G{j)j + cyclic = for a + /3 + 7 = 0, 

which is a sufficient condition for the operator (2) to satisfy the Jacobi identity; the 
general solution to the functional equation was found by Braden and Byatt-Smith 
in the appendix of [i i)]. 

It was observed by Popowicz that, apart from the Hamiltonian Hq, the system 
(4) has additional conserved quantities that can be written as 

Hi = ^ {nm^'^/'^Ym}/^ dx, 

H2 = /(-9n^7i"^m"i/3 + UnxTTixn-^m^^/^ + -4m^m~'^/3) (nTO~2/3^A 

(9) 

where in each case the parameter A is arbitrary. In [Id] it is remarked that, having 
three conserved quantities, the system is likely to be integrable. The existence of a 
mere three (or a few) conservation laws does not guarantee integrability, and a more 
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precise requirement (or better, a definition of integrability) in infinite dimensions 
is tliat an integrablc system sliould liavc infinitely many commuting symmetries 
[13]. In fact, since tfiey contain an arbitrary parameter, eacli of Hi and H2 provide 
infinitely many independent conservation laws for the system. However, a brief cal- 
culation shows that the gradient of each functional appearing in (9) is in the kernel 
of the Hamiltonian operator Z for all A, so that all of these conserved quantities 
are Casimirs for the associated Poisson bracket. Hence, regardless of the choice of 
A, neither Hi nor H2 can generate a non-trivial flow that commutes with the time 
evolution dt- 

The fact that the combination w = nm~'^^^ appears in the conserved functionals 
(9) suggests that it is worthwhile to eliminate either m or n and use this as a de- 
pendent variable. Also, as noted by Popowicz, the conservation laws corresponding 
to Hi are reminiscent of analogous ones for the Camassa-Holm/Degasperis-Procesi 
equations, which provide a reciprocal transformation to an equivalent system with 
different independent variables. We now make use of these observations. 

3. Reciprocal transformation and Painleve analysis. In order to eliminate n 
we can rewrite (4) as 

(mV3), = -(mV3c)„ 
Wt = -Cwx, 

with 

C = 2u + v, m = u~Uxx, wrr?^^ = V — Vxx- (H) 
The first equation is in conservation form, and previous experience with the Degasperis- 
Procesi equation [4] suggests taking the reciprocal transformation 

AX — pda; — Cpdt, -p — m}!^^ 

dT = di, ^^^^ 

so that derivatives transform as = pdx , dt = dr — Cpdx ■ 

In terms of the new independent variables X, T and the dependent variables p, w, 
the system (10) becomes 

ip-')T = Cx, (^3) 
Wt = 0, 

and solving the latter two equations in (11) for u, v we can write 

C = 2u + w = 2m + wm^^^ + {pdxf{2u + w) = 2p^ + wp^ + p{pCx)x- 

Substituting back for Cx from the first of (13) gives C ^ p^ + wp^ + p{p{p~^)t)) x ^ 
and differentiating both sides of the latter with respect to X and substituting for 
Cx once more produces a single equation of third order for p, namely 

PxPxT , (1 -Px)PT , „ , . , r ^ 2 nA\ 

PXXT = \ 5 h 2wppx + [wx + opx)P ■ (14) 

p p'^ 

From the second equation (13), the coefficient w = 'w{X) is an arbitary T-independent 
function of X. It turns out that the presence of this arbitrary function provides 
an obstruction to integrability, from the point of view of the Painleve analysis of 
the partial differential equation (14). It is also easy to calculate the images under 
this reciprocal transformation of the conserved densities corresponding to (9) : the 
density for Hi is transformed to w^, and that for H2 becomes —9'w^~'^w\, both of 
which are trivial (since w is no longer a dynamical variable). 

To analyse the singularities of the equation (14) we apply the Weiss- Tabor- 
Carnevale test. The details of the analysis are almost identical to that for the 
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equation obtained from (1) by an analogous reciprocal transformation, so we will 
only give a brief description of the results and refer the reader to [s] for details 
of a similar calculation. There are two types of local expansion around an arbi- 
trary singular manifold = 4'{X, T) = corresponding to singularities on the right 
hand side of equation (14). For simplicity we can take the Kruskal reduced ansatz 
(j) ^ X — f{T) with / an arbitrary function of T, and it is sufficient to take w to be 
a non-zero constant. For the first type of expansion, p is regular as — > 0, and we 
have p = ±4> + a20^ + as'/'^ -I- . . ., where oii{T) and a'i{T) are arbitrary. The res- 
onances are at —1,1,2, corresponding to the arbitrariness of /, a2, as respectively, 
and all resonance conditions are satisfied, so that this defines the leading part of 
a local expansion that is analytic around = X — f{T) = 0. However, for the 
other type of expansion, we have p ~ X]j>-i Pj't''' which gives resonances —1,2,3. 
The leading coefficient satisfies = — //2, while at the next order /3o = — u'/4. 
The resonance —1 corresponds to /, and the resonance condition for the value 3 
(corresponding to (32) is satisfied automatically, but the resonance 2 corresponding 
to Pi gives the additional condition 

/^-- = -i£-0 (15) 

which is not satisfied unless w = Q (since / is arbitrary). When w = Q the equation 
(14) corresponds to the integrable Degasperis-Procesi equation (see [4]; this is also 
clear from the fact that n = in that case). The failure of the resonance condi- 
tion for non-vanishing w means that the local expansion around a pole must be 
augmented with infinitely many terms in log </>, so that the Painleve property does 
not hold. These logarithmic terms are an indicator that (14) is not integrable, and 
hence the system (4) cannot be. 

4. Travelling waves. Travelling wave solutions of (4) arise by putting 
u{x,t) = U{s), v{x,t) = V{s)^ with s ^ x — ct, 

to get 

M{2U + V -cf ^ Ki, N{2U + V - cf = K2, (16) 
with M ^ U - U", N = V - V", where Ki, K2 are constants, and C(s) = 
2U + V ^ c — fcM~^/^ where k = —K\^'^. From this it is also apparent that 
w — NM^'^/'^ ~ K2K1 ^^"^ ~ I ^constant. Comparing with (12) and (13) is clear 
that travelling waves of (4) are transformed to travelling waves p{X, T) — P{S) of 
(14) moving with speed fc, with 

P(5) = M^/^{s), S - kT, dS = M^l^{s) ds. 

The ordinary differential equation for travelling waves of (14) (with constant w ^ £) 
can be integrated twice to yield 



\ds) ^ k 



p4 + eP'^ + mP^ + cP] +1 = Q{P) 



for k ^ 0, where m is another integration constant. This reduces to an elliptic 
integral of the first kind, 

■ dP 

S + const 
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SO that P is an elliptic function of S. Note that these travelling waves provide 
meromorphic solutions of (14), but this does not contradict the Painlevc analysis 
in the previous section, because for travelling waves the singular manifold is of the 
form (/) = X - f{T) = S - So tor constant Sq, so that f(T) = kT + Sq, implying 
/ = which removes the obstruction to the Painleve property in (15). 
In the original variable s, we have a third kind differential 



pVoiP) 

(so that M(s) = P^{S) has algebraic branch points as a function of s). For partic- 
ular choices of constants, when the quartic Q has a double root, the elliptic integral 
reduces to an elementary one in terms of hyperbolic functions, corresponding to 
smooth solitary wave solutions with the characteristic soliton shape. 

Soliton-type travelling wave solutions must have a constant non-zero background, 
since the requirement that U and V tend to zero as s — > ±cx) implies Ki = = if 2 
in (16), hence fc = and the above analysis does not apply. However, in this case 
we can have a weak solution of (16) which is the peakon solution 

and 

m{x, t) ~ 1a5(x — ct), n{x, t) = 2bS{x — ct), 

where a is an arbitrary constant and c = 2a + & is the wave speed. In the next 
section we extend this to multi-peakon solutions. 



5. Hamiltonian dynamics of peakons. The A^-peakon solutions have the ap- 
pearance of a simple sum of N single peakons but with both the amplitudes and 
positions of the peaks being time-dependent, like so: 

N N 

M(x,t) = ^a,(t)e-l"-«^(*)l, «(x,i) = ^&j(Oe~'""'^^*^', (17) 

where aj{t) and bj{t) are the amplitudes of the waves and qj{t) is the position of 
the peak of both waves. The main result is as follows. 

Theorem 5.1. The Popowicz system (4) admits N -peakon solutions of the form 
( 1 1), where the amplitudes Oj , bj and positions qj satisfy the dynamical system 

N 

dj = 2aj^(2afc-f 6fe)sgn(gj-gfc)e"l«^"«'=l, 

k=l 
N 

bj = bj^{2ak + bk)sgn{qj ~qk)e-^''^-'^''^, 
fe=i 

N 

q, = ^(2afe + 6fc)e-l«^-''^l 
fc=i 

for j = 1, . . . , N . These equations are an Hamiltonian system 
a J = {flj, /i}, hj = {bj, h}, qj = [qj, h} 
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with the Hamiltonian h = 2^j^]^(aj + hj), and the Poisson bracket 



{flj, flfc} = 2ajakSgn{qj - qk)e ™ ' 

= ^sgn(gj - qu) {l - e^l''^ 

{qj,ak} = afce"l«^^*l, 

{aj,hk} = ajhkSgrL{qj - qk)e~^'''~'^'' 




This Poisson bracket has N Casimirs Cj = jor j = 1, . . . , TV. 



The proof of the above resuh, which wiU be presented elsewhere, is based on 
integration of the equations (4) and the brackets (7) against suitable test functions 
with support around each of the peaks. Here it is worth remarking that although 
the phase space has dimension 3A^, fixing the values of the N Casimirs reduces the 
motion onto 2iV-dimensional symplectic leaves. Once this has been done, one can 
eliminate the , say, and solve IN equations for &j , q^ . 

6. Concluding remarks. The evidence from Painleve analysis suggests very strongly 
that the system (4) is not integrable. This raises the question of whether the A''- 
peakon system can be integrable for iV > 1. The Liouville- Arnold theorem requires 
the existence of a further — 1 independent conserved quantities in involution (in 
addition to h and the Casimirs Cj which satisfy {Cj , i^} = for any function F 
on phase space). The first interesting case is the 2-peakon problem, which requires 
just one additional conserved quantity. In fact, a direct calculation shows that the 
independent quantity 



is in involution with the Hamiltonian, {J, /i} = 0, so that the N — 2 peakon system 
is completely integrable. 

There is also the question of whether these peakons are stable solutions. We 
propose to address these issues in future work. 
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